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Abstract 



Let r < PSL(2, C) be a geometrically finite non-elementary discrete subgroup, and let its 
critical exponent S be greater than 1. We use representation theory of PSL(2, C) to prove an 
effective bisector counting theorem for F, which allows counting the number of points of F in 
general expanding regions in PSL(2, C) and provides an explicit error term. We apply this 
theorem to give power savings in the Apollonian circle packing problem and related counting 
problems. 



Contents 



1. Introduction 2 

1.1 Lattice point counting problem 2 

1.2 Sector and bisector count 2 

1.3 Counting F-orbits in cones and hyperboloids 3 

1.4 Acknowledgements 5 

2. Statement of results 5 

2.1 Bisector count 5 

2.2 Application to Apollonian circle packing problem 9 

3. Parametrization 11 

3.1 The group G 11 

3.2 Hyperbolic space 12 

4. Line model 13 

4.1 Complementary series for G 13 

4.2 Lie algebra of G 13 

4.3 Universal Enveloping Algebra oi G 13 

4.4 K-type decomposition 14 

4.5 Useful operators in q 14 

4.6 Normalizing operators 16 

5. Automorphic model 18 

5.1 Spectral decomposition 18 

5.2 Patterson-Sullivan theory. 18 

5.3 Lie algebra elements in KA^K coordinates 19 

6. Decay of matrix coefficients 20 

7. Proof of Theorem 2.1.2 21 

7.1 Setup 21 

7.2 Spectral decomposition 23 

7.3 K-tjpe decomposition 23 

7.4 Main term 24 

7.5 Conjugation invariance 26 

7.6 K invariance 29 

7.7 Bounds in the variables a and a' 30 

8. Proof of Theorems 2.2.1 and 2.2.4 31 

8.1 Setup 31 

8.2 Region of summation 32 

8.3 Smoothing 33 



1 



8.4 Applying Main Theorem 

8.5 Ideal triangle case. . . . 



33 
35 



1. Introduction 



1.1 Lattice point counting problem. Let G be a topological group with a norm || ■ || and 
let r be a discrete subgroup. The question of estimating 



as R grows is known as the lattice point counting problem and was first asked by Gauss for 
G = R^, II • II = II • ||l2, and F = 7?] he showed that 



now called the Gauss circle problem. The current record E{r) — 0{r") for cr = ^ (ignoring 
logarithmic factors) is held by Huxley [17, 18], who improved previous results by van der 
Corput [38] (o" = Iy), Vinogradov (o- = ||), and others. 

Delsarte [4] and Huber [15, 16] initiated the study of this question in the hyperbolic 
setting: G = SL(2, R), F a cocompact lattice, and the norm comes from the distance in H^. 
The difficulty that naturally arises here is that due to hyperbolic expansion: the measure of 
the boundary of the expanding ball grows roughly at the same rate as the measure of the 
ball itself. For this reason it is not clear at the first sight that the main term for N{R) is 
given by the volume of the ball of radius R. Selberg [34] was able to produce an excellent 
error term in the function N{R) for any lattice F using his celebrated trace formula. 

The question of understanding the growth of N{R) for infinite covolume subgroups F 
arose naturally. Patterson [32] and SuUivan [36, 37] developed extremely useful machinery 
for analyzing infinite covolume geometrically finite groups F for G = SL(2,R) and more 
generally for G = SO(n, 1). To each such F they associated a F-invariant probability measure 
on the boundary of the hyperbolic plane (n-space) supported on the limit set of F, and related 
it to the spectrum of the Laplacian on the manifold F\H when the base eigenvalue is at least 
(^^) . Lax and PhiUips [24] used wave equation techniques to produce a very good error 
term for the counting problem in H" for general n and showed that the spectrum of the 
Laplacian on F\H" has but finitely many non-tempered eigenvalues. The main order of 
growth in this case is related to the base eigenvalue S{n — 1 — 5) of the Laplacian on FyH" 
for some 6 = 6r & (^^,n — 1), and equals const ■ e^^ for the number of points in a ball of 
radius R. This is consistent with the fact that the volume of such ball grows hke e^""^)^. 

1.2 Sector and bisector count. A question which arises in applications of lattice point 
counting problem is counting in more general expanding sets, not just in balls. Both in 
the finite and infinite covolume case this problem can be approached by cutting a ball into 

sectors and using them to approximate a more general expanding set. In this setting the 
word bisector refers a subset S* of a semisimple Lie group G = KA^K (known as Cartan or 
polar decomposition) for which there exist contractible Sk-^ G K, Sa C A'^, C K such 



7V(i?) = #{7eF:||7||<i?} 
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that S — Sk^SaSk^- When 5";^^ = K = Sk^-, the set 5" is nothing but a ball. When one 
of Sk^i is all of K. the set S is referred to as a sector. In the case G = SO{n, 1) and 
= K, the set Sk^A'^ C G can be viewed as a subset of H" = G/K, and in the ball model 
this set will be a hyperbolic (and in fact Euchdean) sector, motivating these definitions. 

In the case of lattices in SL(2, R) a very strong result was obtained by Good [9] , and under 
additional assumptions of finite geometric property and critical exponent 5r > | Bourgain, 
Kontorovich, and Sarnak [3] proved a similar counting statement without assuming that F 
is a lattice. A very general asymptotic bisector counting theorem for lattices was proven 
by Gorodnik and Nevo [10]. Oh and Shah [29] proved an asymptotic bisector count for 
G — SO(n, 1) for non-lattices, but their proof relies on measure rigidity and is not readily 
made effective. In the present paper we prove an effective bisector counting theorem for 
G = S0(3, 1) for non-lattices (under additional assumptions). Our approach is most similar 
to that of [3]: the main ingredient is representation theory of L^(r\G). We develop explicit 
i^T-type decomposition for complementary series representations of PSL(2, C) and then use 
the adjoint action of G on conjugates of F to avoid some computation. 

The difficulty that arises in generalizing the current approach to SO(n, 1) is that parts 
of the argument are deeply rooted — albeit less so than that of [3] — in the structure of K- 
types for complementary series representations of S0(3, 1). Many computations use specific 
forms of eigenfunctions, which is the main obstacle preventing direct generalization to higher 
dimensions. Nevertheless we anticipate that an effective bisector estimate should hold for 
SO(n, 1) for all n ^ 4. 

There is another method that can be used to prove bisector theorems (and related limit 
theorems) in the case when F C SO(n, 1) is convex cocompact. In this case the problem 
can translated into the language of symbolic dynamics (cf. [23, 35]); we do not pursue this 
method here. 

1.3 Counting F-orbits in cones and hyperboloids. Let Q be a form defined over Q. The 
problem of counting integer vectors a; G Z" such that Q(a;) = c G Z for a given c is very well 
studied. For definite diagonal forms it is related to Waring's problem [40]. For a general 
definite form one can ask for the growth rate of the number of solutions to Q{x) — c as 
c — > cxD. In the case of an indefinite form one can fix c such that Q{x) — c has a solution, 
and ask for the growth rate of 

#{||x||<x|g(x) = c} 

as X — )■ oo for some norm || • || on Z". When the number of variables n is sufficiently large 
compared to the degree of the form then the problem may be solved by the Hardy- Littlewood 
circle method [39] for certain forms Q. Using different techniques Duke, Rudnick, and Sarnak 
[5] and Eskin and McMullen [7] proved that the number of solutions < X \ Q{x) = c} 

grows like the volume of the corresponding ball (under additional assumptions on Q). 

In a similar fashion one can count solutions to < X \ Q{x) = c} with the 

additional constraint that x lies in a prescribed orbit of a lattice F < Oq(Z), but it is a 
harder question when F is not a lattice. The bisector counting theorem allows us to compute 
H^{x G Txq: \\x\\ < X} with an explicit error term for G = S0(3, 1) and certain groups 
non-lattices F. Note that the varieties {Q = c} are fundamentally different according as c is 
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Figure 1.1: Bounded packing 
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Figure 1.2: Unbounded packing 




Figure 1.3: Apollonian packing in an ideal triangle. 
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positive (one-sheeted hyperboloid) , negative (two-sheeted hyperboloid) , or zero (cone), but 
our approach works in every case. 

A particular instance of this question is the Apollonian circle packing problem. Figures 
1.1 and 1.2 show Apollonian circle packings. For a bounded packing P let N^[T) be the 
number of circles in the packing having curvature at most T. For a periodic packing N^{T) 
is the number of such circles in one period. It has been shown [11, 21] that the N^{T) is 
related to counting points of the orbit of F in the cone on the so-called Descartes form 
which has signature (3, 1). The group F in this case is the Apollonian group (up to finite 
index) which has infinite index in Oq^(Z). We show that 

for constants cp, 6, and e. The power savings depend on the spectral gap for F, which is 
not known numerically. We prove the same counting statement about the packing in Figure 
1.3, which is known as the Apollonian gasket or the packing in an ideal triangle. Moreover 
we can prove effective versions for any circle packing theorem (not just the Apollonian one), 
such as the ones in [30] and [28]. 

1.4 Acknowledgements. I would like to acknowledge the help, inspiration, patience, and 
cheerful optimism of Alex Kontorovich, who introduced me to hyperbolic lattice point count- 
ing problems, guided me in every respect, kindly allowed me to use his pictures, and sup- 
ported my travel through NSF Grant DMS-120937. I am greatly indebted to Yakov G. Sinai 
for acquainting me with dynamical problems in number theory as well as for his kindness 
and support. I greatly appreciate discussions with Peter Sarnak who was always supportive 
and persistent. I would also like thank Stephen Miller for supplying me with vital references. 

I am very grcatful to Ali Altug, Francesco CcIIarosi, and Samuel Ruth for sharing the 
office with me and being always eager to answer my questions. I also thank Arul Shankar, 
Jonathan Luk, Mohammad Farajzadeh Tehrani, Kevin Hughes, Peter Varjii, and other grad- 
uate students who made my time at Princeton mathematically fruitful. 

This manuscript would have been considerably less complete if it had not been scrutinized 
by Hee Oh and Nicolas Templier, whom I thank most kindly. My degree would have never 
been completed without the help of Jill LeClair, who makes simple that which is intricate. 

Finally, I thank my wife, Milena Zhivotovskaya, and my family for their love and support. 

2. Statement of results 

2.1 Bisector count. Let F < PSL(2, C) = G be a non- elementary geometrically finite 
discrete subgroup. The problem of counting the number of points in a F-orbit which lie in 
an expanding region in G is well-studied in the case when F\G has finite volume [34, 9, 5, 7]. 
The main term of such a count is always the volume of the expanding region. There is ample 
literature in the infinite volume case [23, 35, 22, 21, 29]. A major result in this setting is due 
to Lax and Phillips [24], who give an unsurpassed error term using a non-Euclidean wave 
equation. In the present paper we prove a more general counting theorem in the infinite 
volume case. 
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Assume from now on that F is no longer a lattice. The group G naturally acts on 
= {xi + 1x2 + jy I Xi,X2 G R, 1/ > 0} by Mobius transformations (isometries preserving 
the hyperbolic distance d) with quaternion multiplication, and the orbit Tj C is a 
discrete set. The set of its limit points Ap in dM^ is known to be a Cantor set whose 
Hausdorff dimension we denote by 6. As F has infinite covolume we have < 5 < 2 [2]. 
The Patterson- Sullivan measure is a finite measure i^r supported on Ar which is a Hausdorff 
measure of dimension S. This measure is constructed (in essence) as the unique weak-* limit 
of the family of measures 

( \ "^^^ 



■yer 

as s — > 5 from the right [32, 37, 36]. It is the unique finite measure (up to scalar multiples) 
with the property that 



where (3^{zi, Z2) = lim^-j.^ d{zi, z) — d{z2, z) is the Busemann function and 7*^'r = ° is 
the push-forward. 

The group G acts on L'^{T\G): for ^ e G and / e L'^{T\G) we have 

g.f{Tx) = f{Txg), 

which is the right regular representation. Suppose 6 > 1; then the hyperbolic Laplacian 
acts on the smooth functions inside this space and has an eigenfunction with eigenvalue 
Xq — 5{2 — S) < 1. It has but finitely many eigenvalues A„ e (0, 1) with ^ n ^ d, as was 
shown in [24]; arrange the eigenvalues in increasing order. We can write 

for 1 < So < Si ^ . . . < 2 and decreasing in n. As a G-representation the right regular 
representation decomposes as 

where are isomorphic to complementary series representations with parameter s„, and 
So = S. 

In this setting Lax and Phillips [24] prove the following counting theorem. 
Theorem 2.1.1 (Lax-Phillips). Assume S > 1. Then 



1+1 , 



#{7 e F: d{j, 7j) <R} = coe'^^ + Cie^^^ + • • • + Qe^''^ + 0(e— ^i?^). 

The constants Cn are computed explicitly in terms of the corresponding eigenfunctions 
with eigenvalue A„. 

For g E G, let ki{g)a{g)k2{g) be a KA~^K decomposition; it is unique up to factors from 
M, which is the normalizer of A in K. For the group G in question we have 

i^^PSU(2), ^+^R+, M^S0(2). 
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Thus, K/M — and M\K = S"^ with different coordinates. Also define K{g) to be 
ki{g)k2{g). Notice that this definition is independent of the location of the M factor. 

There is an orthonormal basis for L'^{K/M) = L'^{dli^) given by spherical harmonics 
Yabi^, 0), 0, ^ 0, \b\ ^ a; we always assume that every compact group has measure 1. For each 
fixed a the collection {Fab, |^| ^ a,} constitutes a basis for the (2a + l)-dimensional irreducible 
representation of PSU(2). We can, and do, think of spherical harmonics as functions on K; 
in this case they are M- invariant on the right. More generally consider 

L^K)^ J2 (dimT)-T. 

T irrep of K 

Define generalized spherical harmonics 

to be an orthonormal basis for L'^{K) such that for fixed a and c with |c| ^ a 

{Ya-bc- \b\^o] 

is a basis for the left regular representation of K on L'^[K); the usual spherical harmonics 
correspond to c = 0. Writing 



Yaa'-bb'dV, 0, 1p, 6*2, V2) = Ya-bc{V^, 0, 0)Ya'-b'c{^, O2, (^2) ^ 

we have a basis for L'^{K/M.K) given by 

{Yaa'-bb'c{<P:0,'ilj,92,(p2)- \b\: \c\ ^ a; |c| ^ a'}. 

This basis inherits both left and right transformation properties under K. We can write 
Yaa';bb'c{K{9)) using our Coordinates; whenever c = we have 



Yaa';bb'ciKig)) = Yabik^igW^'b' (k^^g)) . 

Finally let | ■ | be a bi-i^-invariant norm on G normalized as explained in (3.1.2) below. 
We prove the following theorem. 

Theorem 2.1.2. Assume 6 > 1. Let 

< 5(2 - 5) = Ao < Ai ^ . . . ^ Ad < 1 

be the eigenvalues below 1 of the hyperbolic Laplacian on r\H^. Then 



+ 



W!>c(i^(7)) = ^^rrr / yab{k)dMk) / Ya,b'{k)diyr{k) 

+ l{e=o}Ci(a, 6, a', 6')^'^^ + ■ ■ ■ + l{c=o}Cd(a, 6, a', 6')^''"+ 

+ 0(r2^(l0gr)^/^i(a + l)i^/^i(a' + 1)1^/^1). (2.1.1) 

The numbers Si, . . . , are less than 5 and > 1 and satisfy s„(2 — s„) = A„. Additionally, 

\crXa,b,a',b')\ < ((a + l)(a' + l))2-^"+^|c„(0,0,0,0)|, (2.1.2) 
and the implied constant depend only on F. 
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Remark 2.1.3. The Patterson- Sullivan [32, 37] measure in the statement is normalized so that 
the corresponding base eigenfunction of the Casimir operator has norm 1; we explain this 
in (5.2.1). 

Remark 2.1.4. The special case with trivial spherical harmonics matches the main term 
Theorem 2.1.1, but the error term in our treatment is worse. 

Remark 2.1.5. The Theorem implies that there are cancellations when c is non-zero. The 
most interesting case, i.e., the case in which we compute the main term, is c = 0. We can 
rewrite it more succinctly as 



ttT r r 

ya'b'{ki{-,'))Yab{k:^\l)) = / Y,,{k)dur{k) / Y^y{k)dur{k) + 

+ ci(a, 6, a', b')T^'' + ■■■ + Cd(a, 6, a', b')T^''^+ 

+ O (T2^(logT)Vii(a + l)i5/^i(a' + 1)1^/") . (2.1.3) 

Remark 2.1.6. The Theorem is only interesting when T is sufficiently large depending on a, 
a'. If the error term dominates the main term a better bound is obtained from Theorem 
2.1.1 by observing that 

\\Ya 'a;b'hc\\L°° < v/(a+ l)(a' + 1), 

whence the right hand side of (2.1.1) can be replaced by 

« T''^/{a + l){a' + l). 



This theorem is a generahzation of the work of Bourgain, Kontorovich, and Sarnak in [3] 
for SL(2, R); we generalize parts of their method to PSL(2, C). In the special case c — 0, the 
main term has been obtained by Oh and Shah [29], but their method uses measure rigidity 
and therefore is not readily made effective. In the present paper we use spectral theory of 
L^(r\G), which provides explicit error terms. For the same reason the method of the present 
paper cannot be applied once 5 ^ 1 as there are no eigenfunctions on r\H'^, but Oh and 
Shah give asymptotics in that case, too. 

The proof of the Main Theorem relies on spectral theory, but we do not follow the origi- 
nal treatment of [3]. We do not attempt to compute the leading term for every combination 
a,b,a',b',c directly — although it is probably possible, and we give tools adequate for com- 
puting the main term in any particular case — but instead develop a new method using 
conjugates of F to obtain the statement for general indices; full computation is only required 
in one case. This new approach works inductively on the indices and relies on the observation 
that once Theorem 2.1.2 is established for one fixed set of indices and all suitable subgroups 
F, the validity of Theorem 2.1.2 can be extended other sets of indices using the freedom in 
F. li g = expeX e G for some X e is close to the identity we can apply Theorem 2.1.2 to 
Fg = g^g~^ and hnearize equation (2.1.1) in e. Equation (2.1.1) for other choices of indices 
can then be recovered as equality between lower order terms for suitably chosen X. 
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A simplified weaker version of Theorem 2.1.2 will be useful to us. Under the assumptions 
of the Theorem, equation (2.1.3) can be replaced by 



^er,|7l<T ^ K/M 

+ O (T'^+%a + l)i5/ii(a' + 1)15/11) . (2.1.4) 

In this version we incorporate the lower order main terms into the error term. 

2.2 Application to Apollonian circle packing problem. Figures 1.1 and 1.2 show integral 
Apollonian circle packings. A thorough introduction to this subject can be found in [11]. For 
a bounded packing P let N^{T) be the number of circles in the packing having curvature at 
most T. For a periodic packing N^(T) is the number of such circles in one period. 

The growth of this function of T was analyzed by Kontorovich and Oh [21] who proved 
the asymptotic formula 

N^{T) ~ cpT\ 

Here cp is some constant depending on P and 5 is the Hausdorff dimension of the packing. 
Since any packing can be sent to any other by a Mobius transformation, the dimension is 
a universal constant S, which McMullen computed to be ^ 1.30568 [27]. The approach of 
[21] is to relate the above problem to counting lattice points in the unit tangent bundle of 
an infinite volume hyperbolic 3-manifold PyTiH^, where F is a finite index subgroup of the 
Apollonian group (see [21, 11]). 
Let be the Descartes form 

QD{a,b, c, d) = + b'^ + + - -{a + b + c + df; 

it has signature (3, 1). To each packing one can associate e with Qd{v) = 0. The 
four entries of v are the curvatures of the four circles (with usual conventions regarding 
signs and orientation) that generate the packing as an orbit of v under the action of F. Let 
V = gu = g{0, 0, 1, 1)"^ for some g G Oq^(R) and let B = g'^Tg; also let B' — {Sig)~^r{Sig) 
where Si is a generator of the Apollonian group (8.1.3). 

In this setting we use Theorem 2.1.2 to prove power savings in the Apollonian circle 
packing problem and compute the overall constant i. 

Theorem 2.2.1. For every e > we have 

r, c 128(5+S1 , 

N^{T) = cpT^ + 0,{T^^+'). (2.2.1) 



^Added in print: Lee and Oh [25] have independently obtained a version of this statement using different 
methods. 
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Remark 2.2.2. The constant cp can be expressed in terms of Patterson-Sullivan measure as 



j (i^|2+i)Wb(^) y 

(B'nAr)\aH3 aH3 

Remark 2.2.3. The constant in front of the main term of this expression is several notation 
changes away from the similar expression in [21]. We use column vectors v instead of row 
vectors, related by transposition. The quantity g from the present paper corresponds to 
{gQ^Y^ in [21], and u is Vq. 

More generally let d: PSL(2, C) — )■ SOq(R) be a fixed isomorphism (used implicitly) for 
a quadratic form Q of signature (3, 1). Let F < SOq(R) be a discrete subgroup subject to 
the hypotheses of Theorem 2.1.2. Let w e {Q = 0} be such that Stabw = NM, let v — gu, 
and let B = g''^^g- Also let || • || be a norm on R^. 

Theorem 2.2.4. For every e > we have 



XG/(stabt,nr)(7) 



||7);||<T 

7er 



6{S-1) J ^11'^ ^ J \\gku 

(BnAr)\aH3 aH3 



128i5+si 

+ 0e{T^2^+'). (2.2.3) 

Remark 2.2.5. In fact || • || only needs to satisfy the scaling property of the norm together 
with non-vanishing and continuity; the full power of the triangle inequality is not needed. 

Remark 2.2.6. We can assume that v e {Q = c} for some c other than zero; a similar result 
holds in this case, but then Stab^; ^ K or Stahv ^ S0°(2, 1). 

Theorem 2.1.2 allows for a further refinement of the Apollonian circle packing problem. 
Consider the packing in an ideal triangle as shown in Figure 1.3. It is contained in the 
region bounded by three tangent circles as shown; it is also contained in the circle through 
the vertices of the triangle. Let 5*4 denote refiection in this circle and let v G R^ consist 
of the curvatures of the four largest circles of the packing ordered so that S/^ preserves the 
three largest circles. Let G4 be the disk bounded by the circumscribed circle of the triangle. 

Theorem 2.2.7. With notation as in Theorem 2.2.1 we have 

r, r 1285+S1 , 

NP{T) = cpT^ + 0,{T^^+') 

where 

^ fn 12 , -,\dj I \ f duB{k) 



G4 9H3 



ii^/c^iiio 



G4 aH3 
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3. Parametrization 



In this section we recall classical facts about the hyperbolic space and related Lie groups. 



3.1 The group G. The group PSL(2, C) is realized as the group of 2 x 2 matrices with 
complex values having determinant 1, modulo the center {±1}. This group is isomorphic to 
SO°(x2+y2 + ^^-«;^) ^ S0°(3, 1). An explicit map from PSL(2, C) to {x"^ + y'^ + -w"^) 
is given by 

a b 



l: ±2\ , I ^ 

/|a|2 + - |6|2 - |c|2 ba + ab-dc- cd_ i{bd -ab-dc + cd) \a\^ + - \c\^ - 

ca — db + ac — bd da + cb + bc + ad i{da — cb + be — ad) cd + db + ac + bd 
i{—cd + db + ac — bd) i{—dd — cb + bc + ad) da — cb — bc + ad i{—cd — db + ac + bd) 
\|a|2 + |c|2 - |6|2 - \d\^ ba + ab + dc + cd i{bd -ab + dc- cd) |ap + |5p + |cp + \d\'^) 

(cf. [11]). We suppress l if no confusion can arise. 

In this realization of S0(3, 1) = SO(a;^ + + — w'^) take the norm 



to be 



Ifl = l«(^)l = 



/cht 



shA 



chy 



It corresponds to 



,t/2 



-t/2 



max(e*/^e-*/^) 



(3.1.1) 



(3.1.2) 



on PSL(2, C). Therefore the statement of Theorem 2.1.2 for S0(3, 1) would have replaced 
by T on the right hand side. 
We set 



e 



-1 



—I 



esu(2,C), h 



1 



-1 



e sl(2, C) 



N = 
M = 



,t/2 



1 z 



-t/2 



: zeC 



-i(f/2 



{exp{th/2): t^Q}, 



■.ipe[0,2n) 



K = {exp((^$/2) exp(^e/2) exp((^2$/2) : G [0, tt), (^2 e [0, 27r), ^ e [0, tt)} . 

The latter parametrization of K is called the Euler angle parametrization. The Haar measure 
on K — PSU(2) in these coordinates is given by 



dk 



47r2 



sin 9 dtp dO d(p2, 
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and the Haar measure on A'^ is 

dat^sh^tdt. (3.1.3) 
On the manifold K/M we get the measure 

dk = — sin 9 dip d6; 

An 

we will use dk to denote the measure both on K/M and K. Thus we have almost everywhere 
unique decomposition 

g = exp((^$/2) ■ exp(^e/2) ■ exp(^$/2) • exp(tV2) ■ exp(^2e/2) • exp((^2*/2). (3.1.4) 

The Haar measure on PSL(2, C) in these coordinates is 

A7rdkidatdk2 (3.1.5) 

Here dki is the probability measure on K/M, and dk2 is the probabihty measure on K. The 
normalization of (3.1.5) is actually quite natural in view of the relation of G to hyperbolic 
space, as we explain below. 

3.2 Hyperbolic space. We use two models for H^, the upper half plane and the ball model 

(see [6] for details). The upper half space model is realized as {xi-\-ix2 + | xi, X2 G R, y > 
0}. The hyperbolic metric and volume are 

2 dx\ + dx2 + dy^ dxi dx2 dy 
ds — 5 dV — 5 . 

y r 

The ball model is realized as the unit ball {xi + 1x2 + jy- + x| + < 1} C R^. We 
use the isomorphism 

{z- j){-]z + l)-^ 

from to to transfer the metric. In these coordinates the metric and volume elements 
are 

^2 ^ '^^1 + '^^2 + '^y^ g dxi dx2 dy 

{1 — xf — xl — y'^y [1 — x\ — X2 — y'^Y 

We also identify boundaries c?H^ and S"^ = K/M by extending the map there. 

In the ball model there is a convenient set of coordinates that match the action of G. 
Any point m'B^ — G/K is the image of under an element of and an element of K/M. 
To wit, 

Wcosf sin|\/e*/2 \ 
e-if/2) [-sinl cosfy 1^ e-*/V"- 

= th -(— cos sin ^ — i sin</?sin 6 + j cos 9). (3.2.1) 

In these coordinates the volume element reads 

dV — sin 9dip d9 ■ sh^ tdt — Andkidat- 

This measure extends to all of G as in (3.1.5). The advantage of these normalizations is 
that respect the structure of G as a X-bundle over H. In particular, the natural map 
L^(H) L'^{G) is an isometry. 
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4. Line model 



The goal of this section is to understand the action of on the i^T-types of a com- 
plementary series representation. Neither the formula from Lemma 4.6.4 nor the action of 
certain operators on the types (4.6.4), (4.6.6), (4.6.5) is widely available in the literature. 
Operators similar to those we introduce appeared in [19] without normalizations. 

4.1 Complementary series for G. Here we review the complementary series representations 
of G in the hne model [20, 8]. Let < s < 2 and let G act on functions /: C — > C by 

Denote by / the intertwining operator 

MC) = Ij^^^^dzAd-z. (4.1.1) 
c 

The action above gives a unitary representation of G on the space of functions from C to C 
equipped with the inner product 

(/i,/2) = J Mz)JM^'-dzAdz; 
c 

call this space Vg. Then the complementary series representation of G with parameter s is 
{tTs, Vs); it is non-tempered and irreducible. 

4.2 Lie algebra of G. We parametrize the (real six-dimensional) Lie algebra g — 5l(2, C) = 
50(3, 1) by 

•=(' -.) "( ') '=(. ) 

and 

*=C _,) 'f-i. )■ 

Note that the letter "i" is concatenated to the symbol it follows. 

4.3 Universal Enveloping Algebra of G. Let denote the complexification of g; it is a 
complex six-dimensional Lie algebra whose field of scalars C we write as a + 61, a, 6 G R. As 
usual, we denote its universal enveloping algebra by U (g^) and lift the representation tt^ to 
it. We shall indulge in the sin of referring to tt^ in this way whether it is a representation of 
G, 0, or U{q^); we shall also drop the subscript s in most cases. 

The center Z of U{q^) is one-dimensional and is generated by the element 

n^h^-ih^ + 2(ef + fe- ieif - ifie); 
it is called the Casimir operator and has the property that 

7r(fi).F + 4s(2-s)F = 

for smooth F & Vs- 
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4.4 K-type decomposition. Consider the decomposition of (tTj, Vs) into K types. That is, 

restrict tt to K < G and decompose Vg into irreducible representations of K. The group 
K has one irreducible representation in every odd dimension. The multiplicity of each 
irreducible representation in 7i\k is 1. Let {vfj-. / ^ 0, |j| ^ /} be an orthonormal basis for 
Vs subject to the conditions 

niQKHj + l{l + l)vt^ ^ (4.4.1) 

7r{ih)vfj = 2ijv^., (4.4.2) 

where Hk = + {ie + ify + {e — fY) is the Casimir element for K. Thus for each I the 
set {vfj : lil ^ 1} is a basis for the {21 + l)-dimensional irreducible representation of K. This 
allows one to recover exact formulae for as solutions to the system of ordinary differential 
equations (4.4.1)-(4.4.2). We will suppress the index s if no confusion can arise. 

4.5 Useful operators in q. The following elements of U{q^) will be important for our 
investigation: 

J± = /±K/ (4.5.1) 

R = ie+if + l{e- f) (4.5.2) 

L^ie + if - l{e- f). (4.5.3) 

The operators R and L (raising and lowering, rcsp.) are the usual ladder operators for 5u(2) 
and map vij to multiples of f/j+i and f^j-i, respectively. 

Lemma 4.5.1. The action of satisfies 

7:{J+)vu = C{l)vi+,,i+, (4.5.4) 

and 

■k{J~)vi^^i = C(/)v,+i 
for certain normalizing constants C{1). 

Proof. We only prove the first statement; the proof of the second one is similar. 

The computation hes entirely in the Lie algebra. Define the following elements of su(2): 

X^C \ Y ^ l\ ^ f I 



-i I V — 1 / \i 

and the following elements from sl{2, C): 

iX= }i iY=( . iZ ^ ~^ 



1) \-i } V-1 

The Casimir operator for su(2) can be written as 

Q^ = l(x2 + y2 + z2). 
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With this normahzation we have 



Tl{VLK)vij = -l{l + l)vij 



for every |j| ^ I. Thus to claim that / + lif is the correct operator we need two conditions 
to hold: 

7r([f]x, / + \^])vu = -2{l + 1)^;^^ (4.5.5) 

and 

7r([X/2, / + \^])vu = !(/ + mvu. (4.5.6) 

This first condition ensures that the Casimir eigenvalue changes from —/(/ + !) to — (/ + 
!)(/ + 2), while the second condition guarantees that the eigenvalue of X/2 is increased by 
I. Since we assume we start with vu, our initial assumptions are 



Tl{VtK)vu = -1(1 + l)Vii, 

n{X)vu = 2llva, 

7r{R)va = 7r{Z + lY)vu = 0. 



(4.5.7) 
(4.5.8) 
(4.5.9) 



Thus we compute 

-^2 + Y^ + Z^ 



-Y -iZ + 1Z -liV 



-{X{iY - liZ) - iX{Y - IZ)) = 



= -{iXl{Z + lY) + {iY - liZ)X + [X, lY - liZ\). (4.5.10) 

The first term contains the raising operator Z + 1Y, so we drop it. The second term contains 
X, which we replace by 211. The last term we massage further to get 

{l + l){zZ + UY). 

Recalling that 7r{Z + lY) = on vu, we can write (4.5.10) as 

confirming (4.5.5). The second condition (4.5.6) is verified directly by observing that 

-Y-iZ + lZ- liY 



X -Y-iZ + lZ-liY 
2"' 2 



= 1- 



□ 



These operators provide an easy way of going from one K representation to another, and 
we call them "jumping operators" for this reason. We have that vij is a multiple of 

7r(L^-^(J+)>oo. 
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In polar coordinates z — re'" we can write the operators as 



and 



For the jumping operators we have 



7r(i?) = e'" [2isr + i{r'^ + \)dr + ~ ~) (4.5.11) 
7r(L) = e"^" \ -1isr - i{r^ + l)dr + ( r - - ) da] • (4.5.12) 



and 



7r(J+) = -e*" ( + -^a ) (4.5.13) 



7r(J-) = -e-'" [dr-^da] . (4.5.14) 



Finally for the last two operators of the basis {J~^, J ,R,L, h, ih} we have 

n{h) = 2s + 2rdr (4.5.15) 

and 

n{ih) = 2da. (4.5.16) 

4.6 Normalizing operators. In this section we establish how the operators J^,R,L defined 
above affect the norm of vij and compute vij. It is well-known that 



t:{R)vi, = 2V(/-j)(^+J + lK.+i (4.6.1) 

and 

n{L)vi, = 2v/(/+j)(/-J + 1)^/^-1- (4.6.2) 

In order to calculate the constants C{1) from (4.5.4) we need to compute vu and the constant 
for the intertwining operator. 

Lemma 4.6.1. We have that vu is a multiple of 

rV'"(l + r2)-^-' 

Proof. Follows from solving the system of differential equations (4.4.1)-(4.4.2). □ 
Lemma 4.6.2. The action of the intertwining operator I : Vg ^ ^2-3 satisfies 



I.vl = (-l)V- 



(r(^^ - 1))^ 



r(/ + s)r(s-/-i)''" 

Proof. It is clear that 



I.Viiij-, a) = const . vf^ *(r, a) 



since / intertwines the action of tIs- The only question is establishing the value of the 
constant. We plug in a = and apply 



dr 

intertwining operator (4.1.1). The result fo^ 



to both sides, recalling the definition of the 

r=0 

lows. □ 
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Lemma 4.6.3. 

Proof. Let 

for 6/ > 0. Then 

C 

and from the previous Lemma we have 

oo 
r-=0 

Evaluating the integral and solving for bi yields the desired result. □ 

The resulting expressions for general vij are given the next 
Lemma 4.6.4. For j ^ we have 

+ s) T{s - I - l)(2l + !)(/- + J-)! 
vij{r, a) = 



Z!7rr(s - 1) 



^ ^ k=0 



A similar formula holds for j < 0. 



Proof. This formula is verified by induction on j. The base case is j — I, which is Lemma 
4.6.3. The induction step is going from j to j — 1; it is carried out by applying the operator 
L from (4.5.12) with normalization from (4.6.2). □ 

Lemma 4.6.4 above allows us to characterize the action of on the X-types. By direct 
substitution we can verify that 
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'-Tr{h)vij = 



7r(J+)t;/j- = 



+ 



7r(J )vij 



+ 



\l + l-s){l-l + s){P - 
(2Z-1)(2/ + 1) 



-vi-i,j - 



l {l + 2-s){l + s){{l + lf- p) 
(2Z + l)(2Z + 3) 



' {I + 2 + + 1 + j){s + l){l + 2 - s) 



(2Z + l)(2Z + 3) 



^{l - - j - -1 + s){l + 1- s) 



(2Z-1)(2/ + 1) 



-'"1-1,3+1 



'{l + 2-j){l + l-j){s + l){l + 2-s) 



(2^ + l)(2Z + 3) 



[21 - 1)(2/ + 1) 



(4.6.4) 



t^z+ij+i - + (4-6.5) 



- V(Z-j + l)(Z+i)7;ij_i+ (4.6.6) 



5. AUTOMORPHIC MODEL 

Now we need to review some representation theory of G in the automorphic model. All 
results are standard, except for formulae from Section 5.3 that are not to the best of our 
knowledge. 

5.1 Spectral decomposition. Consider the decomposition of the right regular representation 
of G on 

L2(r\G) = V^,®V^,®---® V^, © Kemp. 

Each V^^ is isomorphic to a complementary series representation of G with parameter s„ 
and is irreducible; Vtemp is the reducible part consisting of the tempered representations (see 
[8, 20]). Each V^^ contains exactly one positive X-fixed vector, up to scalar multiphcation; 
call this vector These functions are well-defined on r\H^, and satisfy 

Av^n + ■Sn(2 - Sr^Lfn = 0, 

where A is the Laplace-Beltrami operator on the corresponding hyperbolic manifold. 

5.2 Patterson-Sullivan theory. The function (po has the smallest eigenvalue and is called 
the base eigenf unction; it can be realized exphcitly as the integral of the Poisson kernel 
(raised to the power S) against the Patterson-Sullivan measure p. This measure is supported 
on a Cantor set A C dH^ = S"^ and has Hausdorff dimension S (see [32, 37]). That is, the 
connection between u and ipo is 

ifoiif, e,r)= J {P{ip, e, r; v, u)Ydv{v, u), (5.2.1) 
(«,ti)e52 
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where 

1 

Piu^ v V zi] ^ 

' ' ' ' 1 — 2r(sin 9 cos (p sin u cos v + sm9 sin (/? sin u sin + cos 9 cos li) + 

(5.2.2) 

is the Poisson kernel in the Euler angle parametrization. The variable r is the distance in 
the disk model with corresponds to th | in our coordinates on A'^. 

We can also view the hyperbolic three space in the upper half- plane model Xi + 1x2 + jy 
with real Xi, y as in [6]. The action of G on is by Mobius transformations with 
quaternion multiplication. In these coordinates we have 

f f + / ^ 

(Po{xi + 1X2 + jy) = I [— — , , ,2 ) diy{z). 



2ecu{oo} 



\z — (xi + ix2)\'^ + y^ 



There is a similar statement relating eigenfunctions with eigenvalues s„ for n ^ 1 to 
certain distributions -Dr.m not measures, on (9H^ via the Poisson- Helgason transform [12, 1] 
(see [13] for a discussion of the SL(2, R) case). We have 

9, r) = DrAiPiv, 0, r; (5.2.3) 
In [12] it is proven that Dr^n is in the dual space of C^~*"(5'^). 

5.3 Lie algebra elements in KA^K coordinates. We will subsequently need the action 
of U{g^) in KA'^K coordinates; specifically, in coordinates from (3.1.4). This equation 
is a smooth map from to G. Its derivative is a map from Lie(R^) to g. The deriva- 
tive map is invcrtiblc, and we can write the basis {h.ih.e.ie, f.if} in terms of the basis 
{d^, de, d^, dt, de^, d^p^}. We suppress the lengthy computation and only cite the final result. 

^ = — cth t sin ^2^612 + cos ^2^* — ctg 9 csch t sin tp sin 628^+ 

+ cos tp csch t sin 0280 + cosec 6 csch t sin ip sin 92dip (5.3.1) 



— - P) 
2 ~ 



(5.3.2) 



el 1 1 

- = — - (ctg 62 + cth t cosec 62) sin (p2d^2 + 2 (1 + cos 02 cth t)56i2 -|- - cos (p2 sin 025t+ 

1 02 02 

+ -e~* cosec -— csch t sec -r-((— 1 + cos ^2) sm(f2 + e^*(l + cos02) sin (p2+ 
8 2 2 

-|- 2e* ctg 0(cos V' sin (p2 + cos (^2 cos 02 sin ■i/') sin 02)d^+ 
+ ^ csch t(— cos ip2 cos V' cos ^2 + sin ip2 sin ilj)do— 



- cosec csch t (cos ip sin (^2 + cos ip2 cos ^2 sin "0)59? (5.3.3) 
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ie 1 1 

— = — - cos (f2 (ctg 02 + cth t cosec 92)d^^ — - (1 + cos O2 cth t) sin <p2de2 — 

— ^ sin y?2 sin 6281 + ^ (cos (^2 (ctg ^2 cth t + cosec ^2 + cos ■i/' ctg 6' csch t) — 
— cos 62 ctg ^ csch t sin (^2 sin -0)5^ + ^ csch t(cos ?/' cos 62 sin (^2 + cos (^2 sin ilj)dg+ 

+ - cosec csch t(— cos (p2 cos ■i/' + cos 62 sin 932 sin ^)5<^ (5.3.4) 

/ 1 1 1 

— = -(ctg^2 — cth t cosec ^2) sin (^25(^2 + 2 ^'-'^'^2(— 1 + cos ^2 cthf)5e2 + - cos (^2 sin ^2^*+ 

16 9 
+ -e~* cosec — csch t sec —(e^*(—l + cos ^2) sin(/72 + (1 + cos ^2) smip2+ 
8 2 2 

+ 2e* ctg 0(cos sin ip2 + cos (^2 cos 62 sin ^) sin 62)8^+ 
+ ^ csch t(— cos ip2 cos -0 cos ^2 + sin ip2 sin il))de— 

— - cosec ^ csch t(cos 1/^ sin (/?2 + cos ip2 cos 6*2 sin '^l))^^p (5.3.5) 
if 1 1 

— = — - cos ip2 (ctg ^2 — cth t cosec 02)9,^2 + - (— 1 + cos 62 cth t) sin (p2de2 + 

+ ^ sin(y?2 sin^2i?t + cscht(— cos (/?2(ctg^2 + cosec ^2) — 
— 2e* ctg ^(cos (^2 cos ■0 — cos ^2 sin (p2 sin i/;) + e^* cos (p2 tg — )5v~ 
+ - csch i(cos tp cos 62 sin </?2 + cos (p2 sin ip)d0+ 

+ ^ cosec ^ csch t (cos <^2 cos ip — cos ^2 sin (^2 sin -0)^95 (5.3.6) 
These formulae will be useful later. 

6. Decay of matrix coefficients 

In this section wc state well-known results on decay of matrix coefficients. We refer the 
reader to [26] . [41], [21], [14] for proofs. 
Let cu = 1 - e U{q^). 

Lemma 6.1.1. Let (tt, V) be a tempered unitary representation of G. Then for any K-finite 
wi,W2 e V 

{Tr{kiatk2).wi,W2) < te~* dim 'k{K) .Wi\\wi\\l2 ■ •\/dim7r(i^).W2||w2||L2. 
Furthermore, for any wi,W2 & V we have 

{7r{kiatk2).wi,W2) < ie~*||a;('u;i)||L2||(x;(w2)||L2. 
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Lemma 6.1.2. Fix 1 < Sq < 2. Let (tt, y) be a unitary representation of G that does not 
weakly contain any complementary series representation Vs fors > Sq. Then for any K -finite 

wi,w2 e V 

{7r{kiatk2).Wi,W2) < e~^°^^/dim^^{K).Wl\\wl\\L2 • \^dimTr{K).W2\\w2\\L'2- 
Furthermore, for any wi,W2 &V we have 

{'K{kiatk2).Wi,W2) < e~*°*||a;(«;i)||L2||a;(«;2)||L2. 



7. Proof of Theorem 2.1.2 

7.1 Setup. Let 

N{T)= Yl Ya'a-MKil)). 

7er,|7|<T 

Also define frig) on G by 

frig) = Y,,a-MK{g))x\9\<T. (7.1.1) 
Let Ft : T\G x T\G ^ C be defined by 

FT{g,h) = Y,fT{g-'lh). 

Observe that N{T) = FT{e,e). 

For a fixed parameter rj (to be chosen later depending on T) let ijj: G ^ be a smooth 
function supported in a ball of radius rj about e and having integral 1 with respect to the 
measure we fixed in (3.1.5). Let ^ : r\G — > be defined by 

*(x) = 5^V(7^)- 
7er 



Let 



H{T)^ J J FT{g,h)^{h)^{g)dhdg. 



Lemma 7.1.1. 



\H{T)-N{T)\ < (a+l)3/2(a' + l)3/2r2^r7. (7.1.2) 
Proof The left hand side of (7.1.2) is at most 

j j \fT{g-^lh)-fT{i)\^{g)^{h)dhdg. (7.1.3) 

We distinguish three cases. 

1. If I7I > Yz^, then both frig'^lh) and fril) vanish. 
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2. If ^ ^ ItI > then trivially 

IMg-'ih) - /t(7)| « II WfeclU- « + + 

3. Finally if |7| ^ then 

IMg-'lh) - Ml)] <.v\\ya'a;b'bc\\c^ <^ r){{a + l)ia' + . 
Thus the error (7.1.3) is at most 

V(a + l)(a' + l) J] l + 77((a + l)(a' + l))^/' ^ 1 <^ r/Ha + l)ia' + l)f^T 

by Lax-Phillips [24] and Theorem 2.1.1 above. Specifically, we use the leading term of their 
main theorem: 

1 = const •T2^ + 0(T2^-'^), 

7Gr,|7|<T 

where a depends on the spectral gap for L^{T\G). 

□ 

The previous Lemma allows us to work with H rather than with N, and the next Lemma 
massages H into a more pleasant form. 



Lemma 7.1.2. We have 



Proof. We compute 



H{T)= J fT{g){7r{g)^,^)dg. 



FT{g,h)^{h)^{g)dhdg = J2 J J Mg-''yh)^{h)^{g)dhdg. 
ger\G /ier\G '''^'"3er\G h£r\G 

In the inner integral substitute x — g~^^h, whence h — ^gx. Using F-invariance of ^ we get 

Y j J fT{x)^{gx)^{g)dxdg. 

'<^^ger\Gxeg-^'r{T\G) 

Now the result follows by combining the sum over 7 with the integral in x and renaming the 
variables. 

□ 
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7.2 Spectral decomposition. Write ^ = + • • • + + ^temp with each ^„ being the 
orthogonal projection onto Vs„, and ^tcmp is the orthogonal projection onto T^emp- For ^temp 
we use decay of tempered matrix coefficients. We have 

|(7r(A;iatA;2)*temp,*temp)| < *e~*||*temp||^2.2 

from Lemma 6.1.1. The norm W'^'^ is the second order Sobolev norm; it includes the 
norm of the function as well as the norms of all derivatives of orders 1 and 2. Thus we 
get 

J /r(^)(vr(^)*temp, %erap)dg < \\Ya'b'cba\\L4'^ter.p\\w^.2THogT < T^-^^T^logT. (7.2.1) 

Since complementary series representations of G are irreducible, {7T{g)'^sj^,'^g^) vanishes 
whenever Si and S2 are distinct. For the same reason {7i{g)'^si,'^tcmp) = 0. 

To get the simple version of the Main Theorem (2.1.4) we use the decomposition ^' = 
+ where is the projection onto Vs^ and is the projection onto its orthogonal 
complement V^^. In this setting we can use Lemma 6.1.2 to get 



Then we have 



9&G 

and ^'o contributes to the main term. 

7.3 X-type decomposition. It remains to consider the non-tempered part 



J fTig){7iig)^n,'bn)dg. 



" g^G 



Expand ^„ in the orthonormal basis {vij} for the corresponding representation: 

ij 

Here we use vij : T\G — > C to denote a basis in the automorphic model; note that matrix 
coefficients can be computed in any model, while evaluations vij{g), g & G must be done in 
the automorphic model. In KA^K coordinates as in (3.1.4) and (3.1.5) we get 



47r / ^a'a;&'&c(^1^2)X|at|<T(*n,%)('y/'/,*r 



^'j'^''j'K,M\K 
A+,r\G 



Vij{xatk2)vi'ji{xki ) dx datdk2dki. (7.3.1) 
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Now Ya'b'cba a-nd vij have the same ^-transformation properties on the right. Using orthog- 
onahty of matrix coefficients for K we get 



21oKr 



47ry„c(e)i;'c(e)-^^^^-Yy^^^7^PYy j {7r{at)Vac,Va'c)dat. 



t=o 

Finally using the fact that 

{"^n^Vab) = Vab{e) + 0{\VVab{e)\7j) 

we get 

21ogT 

47ry„c(e)l"a'c(e)^^^^jYy(^^7:p^ J {nat)Vac,Va'c)dat+ 

t=0 

21ogT 



+ 0{r]\Vvab{e)\Yac{e)Ya'c{e) ^^^'^^^20/ ^ J {T^{at)vac,Va'c)dat)+ 



t=0 

2iogr 



+ 0{7j\VVa'b'{e)\Yac{e)Ya'c{e) ^ ^J^^^^, ^ J {7r{at)Vac,Va'c)dat)+0{7j). (7.3.2) 

t=0 

The error depends on Lipschitz norms of Vab and Va'b'] we will control their dependence on 
the indices in an indirect way later. 

Remark 7.3.1. A similar expression can be found in Proposition 3.13 of [3], but the signif- 
icance of the Lipschitz norm is overlooked. In Lemma 7.7.2 we show how to control these 
norms, and the same treatment will work in the case of [3] . 

7.4 Main term. In the computation of the main terms of the Theorem we deliberately avoid 
most complications by working out a skeletal form first and then filling in the constants as 
necessary. 

Proposition 7.4.1. We have 



J2 Ya'b'ihi^WabiK'm = CoT'' + CiT'^' + ■■■ + CdT'''+ 



h\<T 



+ O (T2^(logr)Vii(a + l)i5Ai(a' + l)i5/ii j + 0(. . . ) (7.4.1) 

for certain constants Cn = Cn{a,b,a',b'), ^ n ^ d; three dots stand for the error terms of 
equation (7.3.2) with r] as in equation (7.4.5) below. For c ^ we have 

E ya'a-MKii)) = o (r^^(iogr)vn(a + iy^n\a' + ly^f^') . (7.4.2) 

■yer 

|7|<T 
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Lemma 7.4.2. With notation as before the asymptotics for matrix coefficients of {tTsjVs) are 
given by 

(7r(at)<e><'c) = const(e*(^-2)-l^l* + 0{{a^ + a'^ + Ije-^*-!'^'*))- 
Proof. Without loss of generality take c ^ 0. Prom Lemma (4.6.4) we have 

t;,e(r,a) = P2a-c,c(r)(l + r2)-^-<^ 

for a polynomial P2a-c,c- The powers appearing in the polynomial have the same parity; the 
lowest and the highest powers are c and 2a — c, respectively. Combining P2a-c,c(^e*) and 
P2a'-c,c(^) we can write the matrix coefficient as 



Qa+a'-ccl?'^ e*) / r^ + l Y rdr 



r=0 

(5 is a polynomial with indicated maximal and minimal powers of r^, and its coefficients 
depend on e* as inherited from the product -P2o-c,c('"e*) ■-P2a'-c,c('^)- If is natural to substitute 
u — r"^, followed by 

u + 1 

w — 



ue^^ + 1 
These give 



Qa+a' ~c,c (^g2t^]^) 6 ) C^tO 
1 W-e-^'J We^'-J 



2t 



-^ll^ / (7.4.3) 

1 



ignoring absolute constants. Now the integrand is a finite sum of powers of w and can be 
easily evaluated once Q is given. In particular it follows that the main term is of the form 
const • e**""^*~'^* and that the next highest term is 0((a^ + a'^ + l)e~^^~^) from 4.6.4, as 
desired. □ 

Proof of Proposition 7.4-1 ■ It is enough to extract the correct power of T from the integral 

21ogT 
t=0 

from (7.3.2) and to collect the error terms. Using the previous Lemma and the Haar measure 
on A'^ from (3.1.3) we have that 

21ogT 21ogT 

J {n{at)vac,Va'c)shhdt^ j e*^-l^l*cii + 0(r2M). 

t=0 t=0 
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Thus we get 

d 

CnT^'" (7.4.4) 

n=0 

for c different from as in the statement. 

The error term in the Theorem comes from the error contributions (7.1.2), (7.3.2), and 
(7.2.1). 

The optimal choice for t] that makes the errors from (7.1.2) and (7.2.1) equah 

rj = rn(i-'5)(logr)i/^^(a + l)-^(a' + 1)"^. (7.4.5) 

This gives the advertised error term in the statement. □ 

We outhne the calculation of the main term for the case a — b — a' — b' — 0. Plugging 
these values into (7.3.2) gives 

2iogr 

47rt'oo(e)t'oo(e) J {Tr{at)voo,Voo)dat. 
t=o 

We have 

^;oo(e) = i>(0,0). 

For the matrix coefficient we have 

sh(t^ - t) 
(vr(a>oo,t'oo) = ^^33y^, 

which is computed in the line model using Lemma 4.6.3. The leading term (for large t) of 
this expression is 

gt(5-2t 



S-1 

Using the measure from (3.1.3) for the integral we get the leading term of 

TT 



t>(0,0)i>(0,0)T2'^, 



S{S-1) 

as needed. 



7.5 Conjugation invariance. Generalizing this computation seems very computationally 
intensive, so we use a trick to show that once the Theorem holds for the case of trivial 
spherical harmonics we have just considered, it holds in all cases. For g & G, the spectral 
theory of r\G and g~^rg\G is the same. So we consider the discrete group g~^Tg and 
"linearize" in g. In essence wc need to prove a multiplicity one statement for the coefficient 
of the main term by exploiting the transformation property under conjugation. We proceed 
inductively. The base case (constant spherical harmonics) was outlined above. The induction 
step is built on the next proposition. 

Write inv: G ^ G for the inverse map. For functions -Fi, F2, -^3 on K/M, A'^, and M\K, 
respectively, let 

Fi ® F2 ® F,ig) = F^{kM)F2{<9))Fz{k2{g)). 



26 



Proposition 7.5.1. Let Y^a and Ya'a' be a pair spherical harmonics, and suppose by induction 
that ^ 

T™o <^ XT ® Yaa o my{j) = J(JZrY)^^^"'' ' ^^'^^^(0, a). 

Then 

lim 7^ ^ {Ya'a' ® XT ® ^a+l,a+l O inv(7) + Fa'+l,a'+l ® XT ® O inv 

(7)) = 

^^^^ (^'^rla', a')z>r(a + 1, a + 1) + i>r(a' + 1, a' + l)z>r(a, a) j • (7.5.1) 



7er 



Remark 7.5.2. In effect we suppose tfiat tfie main term of Theorem 2.1.2 for one pair of 
spherical harmonics and for every group F, and prove that the main term is right for a new 
combination of pairs of spherical harmonics. 

We need several lemmas to prove this Proposition. 

Lemma 7.5.3. Let i^{z) be the Patterson-Sullivan measure for F and let g & G. Then, the 
Patterson-Sullivan measure for the conjugate of F is 

^g-'rg{z) = P\gj,gz)pr{gz). 

Proof. This is immediate from the properties of the Patterson- Sullivan measure: the measure 
for the conjugated group is a combination of a push-forward with a Jacobian. □ 

Lemma 7.5.4. We have 

P\j+j;^,9)l^^^-Ssm9e'^. 

Proof This is confirmed by direct calculation using (5.3.1)-(5.3.6) apphed to the Poisson 
kernel (5.2.2) raised to the power S. □ 

Proof of Proposition 7.5.1. Approximate the indicator xt in the definition of fx as follows. 
Let be 1 on [0,T], on [T-l-A, 00), and interpolate linearly on [T,T-\-A]. Similarly let 
be 1 on [0, T — A], on [T, 00), and interpolate linearly on [T — A,T]. We will choose 
A < T later. We consider 5+ below; the treatment of 5~ is the same. We have 

Ya'b' ® '^T,A ® ° mv{g'^'jg) = 

■yer 
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-Ug-iTgia, b')Og-irg{a, b)T^' + c^T^^' + • • • + Qr2^'^+ 



5{5-l) 

+ 0(r2T(logr)Vii(a + l)iVii(a' + 1)15/11) + 0{T'-^A). (7.5.2) 

Let > £ > and let g — expsX for some X e g. It is clear that all constants — both 
implied and the constants c„ — are 1 -|- 0{E). 
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Now we expand both sides in e. The left hand side reads 

{Ya'b' <8) 5+^ Yab o inv(7) + 

7€r 

+ eYa'b' ® ^T,A ® ^ab ° inv(7X) + eVa'b' ® O Yab o inv(-X7) + 

+0(£21®S+^®1(7))). 

The two 0{e) terms are similar; we only treat the first one. As t becomes large the only 
derivatives that contribute are 

^ = cos92dt - sin ^29^2 



2 



ell 1 

- = - cos (p2 sin 6281 - - (ctg 62 + cosec 62) sin 992C^v'2 + ^ cos 992(1 + cos 62)002 

ie 1 1 1 

— = -- sin ip2 sin 6*29* - - cos </72(ctg 62 + cosec 6*2)9,^2 - ^^'^ + cos 6*2) sin ip2de2 

/ 1 1 1 

- = - cos 992 sin 92dt + -(ctg ^2 - cosec 62) sin 992f?¥'2 + - cos 932(-l + cos 92)00^ 
z z z z 

tf IX 1 

— = - sin ip2 sin ^2^* - 7: cos </72(ctg ^2 - cosec ^2)9¥'2 + ^(-1 + cos 62) sin (/?29e2 
z z z z 

It is enough to consider the case X = J^. Here we have 

J+ = sin^a^t + i(ctge2 - cosec ^2)6^'^^ 9,^2 + e^'^^-l + 00862)80,. (7.5.3) 

Suppose a = b and a' = b'. The effect of applying J+ is as follows. The first term turns the 
function 

Ya'a' «) o inv(7) (7.5.4) 

into 

T 

-^>"a'a' ® (XT+A - Xt) ® ^a+l,a+l O inv(7), 

where Za+i,a+i{^,(^) = ^aa(<^) ^) sin = const ■ Ya+i^a+i{^,(^)- From (7.3.2) and Lemma 
7.4.1 it follows that 

~2A^ ^a'a' ^ iXT+A - Xt) ^ ^a+l,a+l O inv(7) 

and 

-5 (8) XT (8) Za+i,a+i o inv(7) 

7er 

have identical main terms (when A ^ T), so we use the second form. The other two terms 
of (7.5.3) turn the function (7.5.4) into 

-aYa'a' <8) XT 8) Za+l,a+l O inV. 
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The net contribution from J"*" is 



-(a + d)Ya'a' ® XT «) Za+l,a+l O mV. 

Now we compare to the right hand side. After hnearizing in e we need to study two 
terms: 

J Y,,{ip,e)PSiJ+j;if,e)dMf,d) + J Yaa{-J+M)dyrM. 

The action —J~^{ip,9) is by a Mobius transformation on the boundary Using Lemma 7.5.4 
we rewrite the first term as 



-S J Za+i,a+i{v,0)dur{(p,9). 
Carrying out the calculation gives 



-a J Za+i,a+i{^,0)di^r{<P,0) 
for the second term. This matches the terms on the left-hand side. The term 



Ya'a' ® (g) Yaa O inv(-X7) 

is treated in a similar fashion. 

Since A and e were arbitrary, 0{e) terms on both sides must match, whence the result. □ 

7.6 K invciriance. The previous observation does not allow using completely general com- 
binations of spherical harmonics due to symmetry, so we use K invariance on one side to 
break this symmetry. 

Lemma 7.6.1. Suppose have know that for a finite set I 
A-^^ E E YLj{h{j))Y,,Ak,\j)) = 

— u{L,J)P{l7J^). (7.6.1) 

Then in fact 



7er,|7|<T iL,j,L',j')ex 



^ {L,j,L',j')ei 



7er,|7|<T 
for each (L, J, L', J') G I. 
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Proof. Redefine frig) from equation (7.1.1) by 



Mg) = J2 YLj{hig))xigi<TYL'jik,\g)) 

{L,J,L',J')eX 

and let x & K. Consider the Theorem for the function fT{g>c). We can carry out the proof 
as before until we get to (7.3.1). This equation will now read 



47r J] I YLj{kT^)YL>Ak2)p^j:?A>c)x\aA<T 

T T T I Tl Til 



'J,' J A+,r\G 



(*n, vij) {vi'f, '^n)vij{xatk2)vi'f{xk^ ^) dx dat dk2 dki (7.6.2) 



where p is the corresponding representation of K. That is, each term in the sum transforms 
under K in the same way as 



TT 



J2 HL,JHL',J') 



^ ^ (L,J,L',J')eX 

transforms under K. Repeating this approach on the other side ensures that (7.6.1) holds 
term by term. □ 

It is easy to see from Proposition 7.5.1 and Lemma 7.6.1 that the action of U{q) is 
rich enough to produce every pair of spherical harmonics from the one we checked by hand 
originally, whence the main terms of Theorem 2.1.2 are verified for all pairs of indices. 



7.7 Bounds in the variables a and a'. We confirmed that Co(a, 6, a', b') = ■^^^h'{Ya'b')i^{Yab)- 
By applying the same differential operators R, L, to the Patterson- Sullivan distributions 
(5.2.3) we get that 

Cn(a,b,a',b') = j^^—^Dr,n{Ya'b')Dr,n(Yab). 

Now the bound on the coefficients (2.1.2) in Theorem 2.1.2 follows from the regularity 
property of Dr,n in Section 5.2 since ||Fab||Loo <^ -\/a~+T. 

Remark 7.7.1. A bound in terms of a, a' is also needed in [3], where Theorem 1.5 includes 
the incorrect claim that c„(a, a') are bounded for each n. Nevertheless it can be patched into 
a true statement 

|c„(a,a')| <(|a| + l)'"'"(|a'| + l)'-^". 

For SL(2, R) regularity of automorphic distributions -Dr,n which we use here was analyzed 
by Schmid [33] for all representations and by Otal [31] for the complementary series repre- 
sentations. The results of [3] are little affected by this patch: it can only offset the error 
term slightly for bisector counts involving non-smooth functions. 

Finally we need to establish a bound for Vvab(e) that was needed in Proposition 7.4.1. 
We rely on Theorem 2.1.2, which has been established except for the power of a + 1, a' + 1 
in the error term; currently the error term depends on Vvab{e). The main term has been 
confirmed, so we use it to prove the next 
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Lemma 7.7.2. We have V<b(e) < (a + 1)^-" for s ^ 5 and V^;^^,(e) < (a + l)^-^" for s 5. 

Proof. The action of q on Vab is characterized by equations (4.6.1), (4.6.2), (4.6.4), (4.6.6), 
(4.6.5) above. It follows in particular that for any X e g oi norm 1 we have 

1 1 

i=-lj=-l 

with each Cij -C a + 1. Therefore 

1 1 

7r{X)vab{e) < (a + 1) ^ ^ \va+i,b+j{e)\. (7.7.1) 

i=-ij=-i 

From the fact that the main term of Theorem 2.1.2 is correct we glean that 

Dr,n{Yab) ~ const • Vab(e){a + ly^/^ 
for some absolute constant. The result follows. □ 

This Lemma provides the control needed in Proposition 7.4.1, showing that the two error 
terms there are of the same order in a and a'. The proof of Theorem 2.1.2 in thus complete. 



8. Proof of Theorems 2.2.1 and 2.2.4 

8.1 Setup. Consider the sum 

XG/(Stabi)nr)(7) (8.1.1) 

||7i;||<r 

with V 7^ a column vector in in the cone {Q{v) — 0}. Here Q is a quadratic form of 
signature (3, 1), and assume that an isomorphism i?: PSL(2, C) — >■ SOq(R) has been fixed. 
The norm || ■ || on R'^ can be any norm whatsoever; in fact, it need not even satisfy the 
triangle inequality, but must have the scaling property, some non-vanishing, and continuity. 
For the Apollonian circle packing problem the form of interest is the Descartes form 

Qoia, b,c,d) ^a"^ + b'^ + + d'^ -^{a + b + c + df . 

One can check that 

9 ^ qgq-^ (8.1.2) 

for 





(I 


-1 -1 1\ 


1 


-1 


1-11 


2 


-1 


-1 1 1 






1 1 ij 
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is an isomorphism SO(a;^ + + 



w 



SOq^(R). This map composed with l from 



(3.1.1) defines the action of PSL(2, C) on R"'. Let F be the Apollonian group [21, 11]. It is 
generated by the matrices Si, S2, Ss, S4 given by 



/-I 2 2 2\ 
1 

1 

V 1/ 



/I 

2 

V 



-12 2 
1 



(1 


\ 


1 




2 2 


-1 2 







/I 



(8.1.3) 



^2 2 2 -1/ 



inside Oq^(R). Let F be the image off under the projection Oq^(R) — )■ SOq^(R). Observe 

that these generators have determinant —1, so that products of pairs of generators are in 
SOq^(R). Thus we can count even length words in F directly, and count odd length words 
by looking at the action of F on Siv. In what follows we consider only one of these sums as 
the treatment of the other sum is identical. 

8.2 Region of summation. The set G/ (Stab f flF) should be interpreted as one fundamental 
domain under the action of Stab f fl F on G. It is known that Stab v = NM since v lies in 
the cone {Q — 0}. Since the action of S0(3(R) on the cone is transitive, there are e R^ 
and g E G so that v — gu and Stab u — NM. Then we can rewrite the sum as 

Sv,b{'^) = XG/(StabwnB)(/3) 
||g^w||<T,/3eB 

with B = g'^Tg. Now Stabw = NM, and by an argument from [21] NM n B = iV n B. 
Therefore Stab-uflB is a discrete subgroup of R^ and can be isomorphic to the trivial group, 
Z, or Z^. In the case of Apollonian circle packings only two of these options are possible, 
the trivial group and Z. The first is realized by a bounded packing such as the one in Figure 
1.1; the second corresponds to a periodic packing like the one in Figure 1.2. 

Suppose the set R = Stab -u n B is non-trivial. If G is written in Iwasawa coordinates 
KAN, then R imposes a restriction only in the N variable. That is, R = K.A. Proj^(i?) = 
i^. Proj^^(i?). We denote these projections as R, too. The region is a strip bounded by a 
pair of parallel lines or the intersection of two such strips. In the upper half-space model 
AN = this lifts to a region between two parallel vertical planes or the intersection of 
two such regions. Suppose without loss of generality that j G is within this region. Let 
Ri C H"^ be the region below geodesies joining j to the boundary of ProjjY(-R) C dH^. 
Clearly Ri C R; let R2 be the complement of Ri in R. It is easy to see that the region R2 
is contained within one fundamental domain for F C G. Therefore R2 contains at most one 
point of F, whence 

Il9/3"II<T 
/3eB 

The region Ri conveniently replaces R and is easy to parametrize in KA^K coordinates. It 
is left X-invariant, and in the right K factor it is supported on directions emanating from j 
which hit the boundary inside Proj^(i?). Thus XRiW) = XProig^^s RihiP)) ■ 
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Now in the KA^K decomposition we compute that 



g.kiatk2.u = ^(1 + cos9{k2))\at\.gkiU + 0(1). 

The main term in this expression vanishes (i.e., cos^(/c2) 7^ —1) only if k2U is orthogonal 
to the highest weight vector in the ^-representation. Write 7 = kiatk2 for 7 e F. Then 
cos^(A;2(7)) = —1 for only finitely many 7 e F since its action is properly discontinuous; 
therefore we tacitly assume cos9{k2{'y)) ^ —1. 

8.3 Smoothing. To prepare the sum for an application of Theorem 2.1.2 we smooth the 
indicator xri to approximate from above and from below; let the smooth version Xr^ be non- 
constant in a neighborhood of size U about the boundary of the region S. This introduces 
an error of ^ U^T^. Also take a partition of unity {pY , ^ G /} subordinate to an open cover 
{Vi \ i E 1} with diam(Vi) ^ V for all i G /; we can take |/| ~ V^^. Using an obvious 
approximation argument we recast the sum as 

E E xUk2mprikimpj{k2m. 

''^'^ ^|l+cos9(fc2(^'))lll9fel(/3)«ll 



8.4 Applying Main Theorem. Fix 



W — max /, 

supp unRi 



which exists since the intersection is compact. The function 



^{l + COS0{k2{m 

is unbounded, so we replace it by 

f{k2) A W 

and smooth it at the edge; call the resulting smooth function f{k2). First we consider the 
bounded function / and then estimate the error committed in passing from / to /. 
Pick points Vi &Vi for each i & I. Approximate the above sum by 



E E xlWM{kM)pJ{k2{f3)). 



The error incurred at this step is <^ VT^. 

From Theorem 2.1.2 in the simplified form this sum equals 

Hx''prHpY){m')y ^s ^ o(t^v— 
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The error term comes from estimating Fourier coefficients of x ^-iid P ^"^^ ensuring that the 

sum over a, a' converges. It is clear that the largest error contribution will come from the 
partitions of unity, so we can set [/ = V . Now the sum over i and j approximates the 
respective integral with error <^ VT^ . Thus the leading term becomes 

as in the statement of the Theorem; here we have used the fact that / and / only differ 
outside of the support of v. In the more common coordinates on (?H^ = C U oo the first 
factor reads 

j {\z\' + lfdv^{z). 

(BnW)\(9H3 

The error term comes from optimizing the two different error terms from this section. They 
are 

VT\ T^^T/-^-2(Tf+2)+^. (8.4.1) 

5—s\ 

Equating them gives V — T 129" , which works for all counting theorems stated above. 
Now we estimate the error from bounding /. Sacrificing a constant factor we can focus 

on 



\\gkiil3)u\\ 



■||9fel(/3)«ll 



«E E X{f>w}{k2)pY{h)^Y. E X{f>w}{k2)pl{k2), (8.4.2) 

iel |/3|<r/(fc2(/3)) i m<Tf{k*) 

where k* e 1^ is a point where / attains its maximum over V^. Applying the Main Theorem 
in the simplified form gives 

«E(/((*.')-')7')'^«E''^"'(«v)-'. 

i a 

To ensure convergence in a take I > ^ + 2 and rewrite the expression as 

V^E(/((^n-^))'^^^)-T^^^y- 



.-2-2-if 



The quantity in parentheses asymptotically is a finite integral since / as a simple pole. The 
resulting error term is smaller than both error contributions in (8.4.1), so neglecting this 
term is justified. 
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8.5 Ideal triangle case. Consider the function N^{T) for counting circles in an ideal tri- 
angle, not in the entire packing. Suppose that v G consists of the curvatures of the 
four largest circles in Figure 1.3. A circle from the packing that is contained in a particular 
ideal triangle is characterized by the property that the quadruple in which it appears has 
prescribed form in terms of generators (8.1.3). To wit, the first generator cannot be one that 
flips the smallest circle out of the triangle; say the circles are ordered so that this generator 
is 5*4. The group {1, 5*4} acts on G (and also on K\G = H^) from the right; let 6*4 be the 
fundamental domain for this action that lies under the hemisphere fixed by 154. That is, G4 
contains the ideal triangle in Figure 1.3 by implicitly identifying G4 with its projection onto 
and dH^. The counting function S-^^riT) from (8.1.1) should be changed to 

Sv,AT)- E XG4/{Stab.nr)(7)- (8.5.1) 

||7i;||<T 

Proceeding as in the case of full packings we get the formula of Theorem 2.2.7 with u 
and B deflned as before. 
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